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ABSTRACT 

We present a panchromatic investigation of the partially-embedded, emerging mas¬ 
sive cluster Source 26 (= S26) in NGC 4449 with optical spectra obtained at Apache 
Point Observatory and archival Hubble, Spitzer, and Hersche^ Space Telescope images. 
Eirst identified as a radio continuum source with a thermal component due to ionized 
material, the massive cluster S26 also exhibits optical Wolf-Rayet (WR) emission lines 
that reveal a large evolved massive star population. We find that S26 is host to ~240 
massive stars, of which ~18 are Wolf-Rayet stars; the relative populations are roughly 
consistent with other observed massive star forming clusters and galaxies. We construct 
SEDs over two spatial scales (~100 pc and ~300 pc) that clearly exhibit warm dust and 
polycyclic aromatic hydrocarbon (PAH) emission. The best fit dust and grain models 
reveal that both the intensity of the exciting radiation and PAH grain destruction in¬ 
crease toward the cluster center. Given that the timescale of evacuation is important for 
the future dynamical evolution of the cluster, it is important to determine whether O- 
and WR stars can evacuate the material gradually before supernova do so on a much 
faster timescale. With a minimum age of ~ 3 Myr, it is clear that S26 has not yet 
fully evacuated its natal material, which indicates that unevolved 0-type stars alone 
do not provide sufficient feedback to remove the gas and dust. We hypothesize that 
the feedback of WR stars in this cluster may be necessary for clearing the material 
from the gravitational potential of the cluster. We find S26 is similar to Emission Line 
Clusters observed in the Antennae Galaxies and may be considered a younger analog 
to 30 Doradus in the LMC. 
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Subject headings: galaxies: star clusters: individual: S26 in NGC 4449 — galaxies: star 
formation — HII regions — stars: Wolf-Rayet 


Introduction 


The energetics of galaxies are largely driven by the evolution of massive stars; and when 
clustered, the impact of these stars can be catastrophic. Massive stars modify their environment 
through strong, fast winds and eventual supernova explosions. These pro cesses inject energy an d 
distribute heavy elements, making massive stars critical to galaxy evolution (jMaeder &: Contilll994l '). 
Massive stars are found in the highest concentrations in massive and super star clusters (SSCs), 
which host hundreds to thousands of massive stars in a few parsecs and form in the most intense 
regions of star formation in the universe. However, the physical conditions that produce SSCs, or 
even slightly smaller massive star clusters, remain uncertain, although it is apparent the feedback 
from the constituent massive stars will drive the cluster evolution. In this work, we present a 
massive star cluster undergoing a major transition that may develop our understanding of the 
interplay of massive stars and super star cluster evolution. 

A picture of super star cluster evolution has developed in which SSCs form in thick , dens e 


envelopes of natal material similar to scaled up versions of single massive stars (|Johnsonl 120021 ). 


SSC evolution starts with a molecular cloud proto-cluster that begins to form stars, yet as these 
stars are still embedded in the natal cocoon, the early evolution is effectively obscured from view 
at many wavelengths. However, the massive stars within the embedded cluster begin to ionize 
the surrounding material. A number of these analogs to Ultra-Compact HII (UCHH) regions 


uum sources fe.g. 

Kobulnick'' 

/ &: Johnson 

1999: Turner et al. 

2000|; 

Johnson &: Kobulnickv 2003: 

Johnson et al. 2004 

Tsai et al. 

2ooa 

Reines. Johnson. & CossIboOS: ^ 

Uhnson et al.ll2009: Tsai et al. 

2009: Aversa et al. 

2011: Keolev et al. 201^ 

) with a ffat or inverted spectral index, indicative of 


thermal free-free emission from dense young HII regions. iKobulnickv fc JohnsonI (119991') dubbed 


these sources as Ultra-Dense HII regions (UDHHs) (or similarly “supernebulae” by iTurner et al 

2nnnl b 


The massive stars will continue to evolve and proceed to evacuate the surrounding material. 
One possible example of this emerging evolutionary stage are emission line clusters (ETC). Iden¬ 
tified in the Antennae galaxies, ELCs are a type of HII region that are younger versions of SSCs 
and exhibit broadened Bry line emission suggestive of massive stars evacuating their surroundings 


^Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator 
consortia and with important participation from NASA. 
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via wind (jGilbert &: GrahamI l2007l i . Finally, the massive stars will be revealed at optical wave¬ 
lengths and regulate or halt star formation (jAgertz et alJ [2013) • This results in the final early 
evolutionary stage of SSGs a s bright and blue optical clusters that are well studied with HST 
([Whitmore &: Schweizen Il995l ). This last stage can be exemplified by the well-known region 30 
Doradus (30 Dor) in the LMG, which is the closest SSG analog. The entire massive/super star 
cluster evolutionary sequence can be summ arized as: protocluster —)• UDHII —emerging cluster 
(ELG) —7> SSG (e.g. Whitmore et ah 2014 ). 


Although some of the steps in the evolution have been outlined, ultimately the physical pro¬ 
cess of an UDHII region becoming a cleared-out, optical SSG is not yet well understood. For 
instance, cluster age measurements are not feasible until the stars are optically visible after clear¬ 
ing embedding material and the most massive stars have started to evolve off the main sequence. 
Popu lation comparisons suggest that the UDHII phase lasts for < 1 Myr ([Kobulnickv fc Johnson 
19991 1 ■ however, a radio and o ptical study of NGG 4449 reve aled that some clusters may in fact re¬ 
main embedded up to 5 Myr ([Heines. Johnson, fc GosaDOOSi ). As for the ev acuation process, many 


mode ls assume an instantaneous removal at some given age (such as in IPfalzner fc Kaczmarek 


2 OI 3 II. This picture is too sim plified, as the rate of removal will surely change the fate of the cluster 
(iBaumgardt &: Kroupall2007ll . E ven accounting for the rate of removal is not enough: as shown 
by jPfalzner &: Kaczmarek! ([20131 ) , a scenario simply comparing SEE or expulsion timescales is too 
limited to describe a cluster’s ability to stay bound and thus survive. 


Perhaps most paramount, the dominant mechanism responsible for the evacuation of the natal 
material is unclear. Massive stars erode the obscuring envelope through a combination of stellar 
feedback processes including direct radiation from stars; pressure from cold, warm (the ionized HII 
region itself), and hot gas; dust processed IR radiation; protostellar winds and jets; and stellar 
winds and supernovae (e.g. iLopez et al.ll2013l ). Yet, the relative importance of these mechanisms 
in removing the natal material is under d ebate - especially how t hese processes are coupled to the 
molecular cloud material (as discussed in [Rogers fc Pittardll2013lb 


We have identified a young massive cluster in NGG 4449 that appears to be undergoing this 
major transition between evolutionary phases, and its stellar content is consistent with winds 



drastically impact their environment. Thus, we hypothesize that we may have caught S26 in the act 
of breaking out of its natal cocoon, driven by the winds of the WR stars, providing an opportunity 
to observe the impact of evolved massive stars on their natal environment. 


S26 is in a relatively low metallicity environment, and thus not only represents an important 
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Fig. 1.— A Hubble Space Telescope rgb image (Ha, I, B) showing nebular emission surrounding 
a compact optical cluster in S 26. The green contours show 3,4,5, and 6a emission at 3.6 cm 
(jReines. Johnson. &: GossI 1200811 : the 3 cj contour corresponds to a region with a radius of ~50 pc. 
Color version available online. 


stage in cluster evolution, but also offers a window into a regime of star formation that is not well 



We present a multi-wavelength analysis of the massive star cluster S26 in NGC 4449. The 
optical spectra and archival infrared data and reduction are discussed in Section [2j We evaluate 
general properties of the region, such as extinction and age, in Section [3l We determine the massive 
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Wavelength [A] 


Fig. 2.— Flux calibrated spectra taken with DIS on the 3.5 m Telescope at Apache Point Obser¬ 
vatory of the cluster S26 in NGC 4449, zoomed in on the broad Wolf-Rayet features. The dotted 
line shows the subtracted nebular features. Top- The “blue bump”, a composite of broad lines at 
4640 A, 4650 A, and 4686 A. Bottom- The broad “red bump” centered at 5808 A. Obvious nebular 
and WR features are labeled. 


star populations in Section [4] and identify the thermal radio emission in Section [5l In Section [6l 
we construct SEDs and find the best fit dust models. In Section [71 we put S26’s massive star 
populations in context, discuss winds as possibly driving the evolution, and look at the similarities 
between S26 and 30 Doradus in the LMC. Finally, we present our conclusions briefly in Section [8l 


Observations and Data Reduction 


2.1. Optical Spectra from APO 


S26 in NGC 4449 was obs erved on 2008 April 13 with the 3.5 m telescope at Apache Point 
Observatory (jReines et al.ll2010l ). using the red and blue channels of the Dual Imaging Spectrograph 
(DIS) in low-resolution mode. The total exposure time was 30 minutes, producing a signal-to-noise 
ratio per pixel of ~60 in the blue and ~45 in the red continua. The spectrum has a resolution of 
~7 A over a wavelength range of ~3800-9800 A and was reduced using IRAF and IDL routines. 
S26 was observed with a 1”.5 x 360” slit and its spectrum was extracted from a 4”.4 window. The 
correction factor to account for slit loss is 1.9, determined by co mparisons to HST p hotometry of 
a region of radius 3”.3 (~ 60 pc). Further details can be found in iReines et al.l (|2010j). 


The spectrum of th e cluster S26 clearly displays typical Wolf-Rayet line features (Fig. |2l 


also iReines et al.l 120101 1. explained in detail in Section [41 After the initial spectral reduction 
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Fig. 3.— An infrared view of the massive cluster S26 in NGC 4449. Prom top left moving clockwise: 
an Spitzer IRAC rgb image (3.6 pm, 5.8 ^m, and 8.0 ^m), Spitzer MIPS 24 pm, Herschel PACS R 
(170 pm), and Herschel SPIRE PLW (500 pm). The extraction regions for constructing the source 
SED (Eig. E]) are plotted in red (small 5”~100 pc circle) and green (large 15”~300 pc circle). S26 
is a dominate source of emission even at the longest wavelengths. Color version available online. 


of iReines et al.l (j2010l l. emission line features are further processed in this work using the IRAE 
SPLOT package. The blue bump WR feature, a composite feature of lines at 4650 and 4686 A, 
is analyzed after the prominent, superimposed nebular lines [Ee III] (4658 A) and Hel and [Ar 
IV] (4713 A) are subtracted (Eig. [2]). The observed line flux, extinction corrected line fluxes (see 
Section 13.ID. and equ i valent widths are given for emission lines in Table [TJ Results are similar to 


those in 


Reines et al.l (|2010l l for the only presented line, Ha. 













Table 1. 


Emission Line Properties of S26 


Wavelength 

(A) 

Identification 

Observed Flux 
ergs cm~^ 

Extinction Corrected Flux 
ergs cm~^ 

EW 

(A) 

3835 


0.195 (0.013) 

0.306 (0.059) 

7.8 (0.7) 

3868 

[Nelll] 

0.754 (0.033) 

1.18 (0.22) 

29.5 (2.7) 

3889 

HC 

0.527 (0.023) 

0.82 (0.15) 

21.1 (1.5) 

3970 

He 

0.607 (0.026) 

0.95 (0.17) 

25.1 (2.0) 

4076 

[SII] 

0.019 (0.013) 

0.030 (0.021) 

0.8 (0.6) 

4102 


0.720 (0.023) 

1.12 (0.19) 

33.5 (1.5) 

4341 

H7 

1.323 (0.042) 

2.03 (0.33) 

70.0 (4.1) 

4363 

[OIII] 

0.066 (0.017) 

0.101 (0.031) 

3.6 (1.2) 

4471 

Hel 

0.116 (0.011) 

0.177 (0.033) 

7.1 (1.0) 

4650 

[CIII] (blue bump) 

0.102 (0.018) 

0.155 (0.037) 

7.1 hs)"" 

4658 

[Felll] 

0.015 (0.001) 

0.0223 (0.0035) 

0.94 (0.02) 

4686 

Hell (blue bump) 

0.077 (0.013) 

0.116 (0.026) 

5.3 (0.9)=" 

4711 

[ArIV]/He I 

0.014 (0.001) 

0.0213 (0.0033) 

0.96 (0.02) 

4861 

H/3 

3.00 (0.11) 

4.52 (0.68) 

185 (49) 

4959 

[OIII] 

4.33 (0.24) 

6.50 hoO) 

300 (200) 

5007 

[OIII] 

13.40 (0.46) 

20.07 (2.96) 

940 (750) 

5755 

[Nil] 

0.009 (0.001) 

0.0125 (0.0017) 

1.0 (0.1) 

5808 

[CIV] (red bump) 

0.040 (0.026) 

0.059 (0.039) 

4.9 (3.4p 

5876 

Hel 

0.348 (0.013) 

0.510 (0.067) 

43.6 (3.6) 

6300 

[OI] 

0.0410 (0.0035) 

0.0595 (0.0087) 

6.13 (0.72) 

6312 

[SHI] 

0.0486 (0.0033) 

0.0705 (0.0096) 

7.25 (0.70) 

6548 

[Nil] 

0.101 (0.070) 

0.15 (0.10) 

8.7 (8.5) 

6563 

H« 

9.60 (0.32) 

13.86 (1.66) 

580 (380) 

6584 

[Nil] 

0.350 (0.093) 

0.51 (0.15) 

46 (28) 

6678 

Hel 

0.106 (0.0051) 

0.153 (0.019) 

18.9 (1.6) 

6717 

[SII] 

0.294 (0.024) 

0.423 (0.059) 

59 (27) 

6732 

[SII] 

0.218 (0.023) 

0.314 (0.049) 

44 (23) 

7065 

Hel 

0.0778 (0.0066) 

0.112 (0.015) 

15.9 (2.7) 

7137 

[Arlll] 

0.344 (0.013) 

0.493 (0.056) 

66.5 (8.6) 

7319 

[OH] 

0.0728 (0.0039) 

0.104 (0.012) 

16.1 (1.5) 

7330 

[OH] 

0.0584 (0.0034) 

0.083 (0.010) 

13.0 hs) 

7751 

[Arlll] 

0.0827 (0.0048) 

0.117 (0.014) 

19.3 (2.3) 

9069 

[SHI] 

0.974 (0.043) 

1.36 (0.14) 

500 (630) 

9532 

[SHI] 

2.503 (0.086) 

3.49 (0.33) 

b 


Note. — Tabulated quantities are: intrinsic wavelength, common identification, observed fiux, extinction 
corrected flux, and equivalent width (EW). Uncertainties follow in parentheses. 

^The Wolf-Rayet features are broad. 

^Not well constrained. 
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Infrared Archival Data 


2.2.1. Spitzer IRAC and MIPS 


Infrared data from the Spitzer Space Telescope (|Werner et al.ll2004l i were retrieved from the 
Spitzer Science Center Archive, consisting of InfraRed A rray Camera (IRA C) (jFazio et alJ 12004 ) 
and Multiband Imaging Photometer for SIRTF (MIPS) ( Rieke et al. 2004) imaging. NGC 4449 
was imaged with IRAC (PI: G. Fazio) at 3.6 /rm, 4.5 /rm, 5.8 /rm, and 8.0 um (FWHM of ^1.9”; 


Fazio et al.l l2004l i and MIPS (PI: R. Kennicutt) at 24 /rm (FWHM of ~6”; iRieke et al.l l2004li 


We analyze the post-basic calibrated data from IRAC and MIPS images, w hich were reduced w ith 
the Spitzer Science Cente r pipeline. The flux cali bration uncertainty is 2% ([Reach et al.ll2005l i for 
IRAC and 4% for MIPS ( Engel hr acht et al. 2007 1. 


2.2.2. Herschel PACS and SPIRE 


The ESA Herschel Space Observatory 

(Pilbratt et al. 

2010') has observed NGC 4449 as part of 

the Dwarf Galaxy Survey ( 

Madden et al. 

2 OI 3 I). NGC 4449 was observed with the Photodetector 

Array Camera and Spectrometer (PACS) 

(Poglitsch et al. 

20101 at 70. 100. and 160 pm ('FWHM 


of the PS F is 5.2. 7.7. and 12 .0”. respectively) and the Spectral and Photometric Imaging REceiver 
(SPIRE) ( Griffin et ahlboiO lat 250, 350, and 500 /rm (FWHM of 18.2, 24.9, and 36.3”, respectively). 


We analyze the newly released ‘Level 2.5 ’ data from PACS, using MadMap images to preserve 
extended stru cture (IRemv-Ruver et al.ll2013l ). The PACS calibration uncertainty is of the order 
of 10%, e.g. Fritz et al. (2012). For the SPIRE data, we analyze the ‘Level 2.0’ and adopt a 
conservative uncertainty of 15% on flux calibration, as according to the SPIRE observer’s manual 
(http://herschel.esac.esa.int/Docs/SPIRE/html/) the overall calibration uncertainty for the 
SPIRE photometer is 7-15%. For both PACS and SPIRE datasets, we apply astrometric corrections 
to the data using the peak emission of the nucleus of NGC 4449. 


3. General Properties of S26 
3.1. Extinction 


The extinction of a star forming region is crucial for accurate line measurements, as well as in¬ 
formation on the extent to which a source is embedded. Using nebular lines only measures extinction 
towards gas that is not very heavily extincted and thus can be biased lo w. We derive the extinc¬ 
tion o f S26 through optical nebular Balmer line and radio fluxes fr om iReines. Johnson, fc: Goss 
(120081). Extinction cu rves for the 30 Doradus region of the LMC (jMisselt. Clavton. fc Gordon 


1999 


Fitzpatricklll985l ) are used to convert the Balmer decrement to Ay, appro priate as the dust 


processing and the metallicity of the two regions are similar, roughly z=0.008 ([Russell fc Dopita 
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1990|)- Use of a Milky Way extinction curve changes the measured extinction by only ~10%. 


The total extinction of S26 is measured to be Ay = 0.41 and is used in correcting the optical 
line fluxes in Table 1. Subtraction of the Galactic foreground extinction results in an internal 
extinction of Ay, i = 0.35: along the S26 line of sight, the Milky Way galactic extinction is measured 
as E(B—V)= 0.019 (ISchlegel. Finkbeiner. &: Davislll998l ). which we convert assuming the standard 
galactic curve Ay = 3.1 E(B—V). Within the un c ertain ties, this internal extinction is in agreement 
with Ay^ i = 0.40 as estimated by iReines et al.l (j2010l l by fitting the spectrum with Starburst99 
models. 


3.2. Electron Temperatures, Density, and Pressure 


In order to probe the pressure as well as measure oxygen abundance as a proxy for metallic- 
ity, we estimate the electron density and temperatures in S26. These physical conditions in the 
i onized gas of S26 are d etermined through ratios of line fluxes using the five-level atom model 


(jPe Robertis et al.lll987l ) with the NEBULAR package in IRAF. The electron density is estimated 


using the SII line ratio 6716A/6731A and the S"*" electron temperature given by the line ratio 
(6716A-|-6731A)/4076A. The electron density is estimated to be ne(SII)Rs75 cm“^, although it may 
vary by as much as a factor of four due to propagation of flux uncertainties also impacting the tem- 



tio (6548A+6584A)/5755A. The electron temperature is measured to be T(OIII) = 9400 ± 
500 K using the [OIII] ratio (4959A+5007A)/4363A. Using the estimated density of 75 cm“^ then 
implies a pressure of P/k = 7.5 xlO^ cm“^ K for S26. The temperatures, density, and pressure 
estimated for S26 fit t he observed range of typical HII regions in disk g alaxies and ELC regions in 


the Antennae Galaxy (iGilbert fc GrahamI 120071 : iHunt fc Hirashitall2009l ). We assume the standard 


HII temperature of 10“^ K as a representative single temperature throughout this work. 


3.3. Metallicity via Oxygen Abundance 


We derive the oxygen abundances using the s tandard Tp method with two distinct temper¬ 
ature zones in the photoionized HII region, as in lizotov et al.l (|l994l . 119971 ). using the O"*" and 
0“*"^ electron temperatures are explained above. The total oxygen abundance is derived by 0/H 
= 0+/H++0++/H+: 0+ ionic abundance is measured with the summed flux of the doublet 
7319/7330A, and ionic abundance is determined using ionic abundances for the lines 4363A, 
4959A, and 5007A. We measure a total oxygen abundanc e of 12 -|- logfO/H ) = 8.3 ± 0.2, which 
is in agreement with checks using empirical relations from lizotov et al.l (j2006l i that were computed 
with new photoionization models. By assuming a simple scaling relation and the solar metallicity 
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Fig. 4.— Starburst99 predictions for the cluster age versus the equivalent width of the H/3 (4861 
A) line (see Section [7.1.11) . The observed EW of H/3 for the massive star cluster S26 in NGC 4449 
(solid line with uncertainties shown as a shadowed region) implies an age of 3.1 ± 0.3 Myr. Color 
version available online. 


value 12+log(O/H)0 = 8.69 (jAsphmd et al.l l2009l ) . the oxygen abundance converts roughly to a 
metallicity of Z~0.006 (=0.4 Zq), similar to th e LMC. The m e asure d oxygen abundance of S26 is 
in excellent agreement with the value found bv iLeaueux et al.l (|l979li of 12+log(0/H) = 8.38. 


3.4. Age 


Ages are instrumental to characterizing the timescales of the early evolutionary stages of 
massive star clusters, as well as necessary for comparisons to predictions and other observed objects. 
The age of S26 is estim ated by compar ing the equivalent width of H/3 to the evolutionary synthesis 
models of Starburst99 ( Leithereiill999l i. using the measured equivalent width of H/3 of «185 A, and 
adopting a metallicity of Z=0.008 (see details provided in SectionlT.l.lD. A s shown in Figure [H we 


derive an age of 3.1 ± 0.3 Myr for S26, consistent with lReines et al 


( 120101 ). 


Wolf-Rayet Features and the Massive Star Populations 


The first hint of the import ant evolutionary p hase of S26 was the discovery of an optical 
feature due to Wolf-Rayet stars (jReines et al.ll201Cl b WRs are spectroscopically identifiable via 
unique emission lines and are divided into subtypes of nitrogen-rich WNs, carbon-rich WCs, and 
oxygen-rich WOs; the presence of several WR stars will produce broad features in the integrated 
optical spectrum known as WR “bumps” (Fig. [2]). S26 clearly displays broad WR features of both 
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WN and WC populations, namely the “blue bump” near 4650 A, consisting of a blend of [NIII] 4640 
A, [CIII] 4650 A, and Hell at 4686 A, and the “red bump” due to [CIV] at 5808 A. The emission 
at 4686 A of S26 does not display a clearl y dist i nct ne bular component (Fig. [2] ) often seen in 
the integrated spectra of WR galaxies e.g. IContil (jl99lh . However, at relatively lo w metallicities, 
such as S26, the stellar component is thought to dominate (jSchaerer &: Vaccalll998l l. The expected 
ne bular emission for a popu lation of 3 - 3.5 Myr should be log( hnebuiar^Pfen [4686 A]) according 


to 


Schaerer fc Vaccal (Il998l ). Thus, the expected nebular emission is too low to be signihcantly 


detected or distinguished from the stellar emission at 4686 A with this dataset, and thus the lack 
of resolved nebular emission is not surprising. Other WR features that could further classify the 
WR population would be similarly too faint to detect. 


4.1. Determining the Number of WR stars 


The luminosity of 
mating the WR star poi 

a WR feature can be used to constrain the WR population and esti- 
pulations in an extragalactic source has become fairly standardized (e.g. 

Schaerer & Vacca 

1998: 

Guseva. Izotov. &: Thuan 

200 c 1. The number of WR stars can be esti- 


mated as Vwr = TwR/-bo,WR) with To,WR as a typical single WR star producing the feature and 
Lwr as the observed luminosity of that feature from the source. The subtypes can be separately 
analyzed with their respective WR lines. Many uncertainties result from this method due to the 
inherent range of WR line fluxes emitted by individual stars, averaging over WR subtypes, and the 
currently poor understanding of the impact of environment on WR line fluxes. Nonetheless, it is 
useful to estimate the massive star populations within a cluster, as long as the range of uncertainties 
are kept in mind while interpreting the results. 


In order to constrain the WR subtype populations in S26, we adopt the method of ICuseva. Izotov, fc Thuan 


( 2 OO 0 I ) . The red bump at 5808 A is thought to be produced only by WCE stars and can be com¬ 
pared to representative WC4 stars. The number of WC stars can be approximated as V-wc = 
L(5808 A)/Lwc4(5808 A), assumi ng a typical luminosity of a single WC4 star at 5808 A is 3.0x 10^® 
ergs s“^ as measured in the LMC (jSchaerer fc Vaccal Il998l ). 


Determining the number of WN stars is less straight forward because the blue bump includes 
contributions from both WN and WC stars. The relative contribution of WC stars to the blue 


coefficient k = Lwc4(4650 A)/Lwc4(5808 A) ( 

Guseva. Izotov. V Thuan 

20001. We adopt a value 

of k=1.71 ± 0.53 ( 

Schaerer V Vacca 

1998: 

Guseva 

Izotov, h Thuan 

200f)ll. although uncertainties 

are large due to variations in relative line fluxes ( 

Schaerer V Vacca 

1998|l. After subtracting the 


estimated WC contribution from the total measured flux in the blue bump, the number of WN 
stars is found by comparing the remaining emission to a typical WN star. We assume a typical 
WN of a WNL (WN7) star to be with a lu minosity of 2.0 x 10^® erg s“^ in the blue bump (4650 
+ 4686 A; Guseva. Izotov. &: Thuan 2000 1 as in the LMC. 















































- 12 - 


In total, we have estimated 20 ± 14 Wolf-Rayet stars in S26, comprised of roughly 4 ± 3 WC 
and 16 ± 13 WN stars. The results can be found in Table [2j The uncertainties are estimated from 
assum ing a typical flux, where we use the observed flux range in the LMC by ISchaerer &: Vacca 
(| 19981 1 ■ as well as from flux measurement uncertainties in Tabled) Due to the large range of observed 
variation of fluxes, the uncertainties are correspondingly, and unavoidably, large. Another potential 
uncertainty that has not been accounted for arises from the weak Hell line emission from the most 
massive 0-stars (Of-stars), which may result in an overestimate of WN stars. However, this effect 
is expected to be small at the metallicity and age of S26, roughly Of/O~0.15, w hich with the weak 
Hell emission (10% of a WN star) would contribute a total of ~2 WN stars (jSchaerer &: Vacca 


199a). Additionally the ratio of WC to WN stars appears normal, see Section 17.1.21 and thus 


Of-star contamination is likely not important. 


4.2. Determining the Number of O-stars 


Identifying the total massive star population in S26 is important for understanding the feedback 
processes that are altering the cluster. In addition to the massive WR stars that are driving 
strong winds and hard radiation throughout the cluster, the 0-star populations will also contribute 
ionizing radiation and additional (albeit weaker) winds. The number of O-stars can be estimated by 
determining the population needed to produce the observed ionizing photons Qo, after subtracting 
off the Wolf-Rayet contributions. The ionizing flux s een at optical wavelengths is simply estimated 
through empirical relations from ISchaerer &: Vaccal (199a). Using H/3 at 4861 A, the ionizing flux 
is Qo~170xl0^® photons s“^. We can thus estimate the number of O-stars by assuming 


= (Qo - AwrQo,Wr)/(??oQo,07v) 


(jOuseva. Izotov, fc ThuanI 120061 ). We include the parameter rjo as the ratio of 07V stars to all O- 
stars to account for different 0-star subtypes occurring within IMF (as subtypes produce different 
i onizing photon fluxes). At an age of 3.1 Myr and Salpeter IMF, we find r]o^l.2 for the cluster S26 
([Schaerer &: Vaccal 1 1998l i. near the peak value resulting from a WR rich phase. A typical ionizing 


Table 2. Massive Star Populations in S26 


Type 

Number 

O 

219 ± 15 

WNL 

15 ± 12 

WCE 

3 ± 3 

WRtotal 

18 ± 13 

WR/O 

0.084 ± 0.058 

WC/WN 

0.23 ± 0.26 
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s ^ from an 07V star (Martins et al. 


2005|) 


photon flux for an 0-star is Qn.mv is taken to be 10 __ 

and 10^® s“^ for a WR star ( Guseva. Izotov, fc ThuanIboOO : SAaerer. Contini. fc PindaJl999 h We 


estimate 0-star population uncertainties by accounting for the measured flux uncertainty in S26, 
estimated uncertainties on ijo, and uncertainties for the subtracted WR populations. With a WR 
population of A^wR = 18 from above, we find there are approximately 219 ± 15 total 0-stars 
harbored in the massive star cluster S26. 


5. The Thermal Radio Component of S26 


The radio spectral indices of S26 indicate that the source hosts mixed thermal and non-thermal 
contributions; and as shown in Figure [H the radio continuum contours additionally change from 
a resolved, irregular shape at 3a to an unresolved, compact peak at 6a. Despite the messy na¬ 
ture of S26, it is necessary to identify the radio emission produced by the HII region itself in S26. 
We d ecompose the observed radio em ission to the published fluxes at 1.3 cm, 3.6 cm, and 6.0 cm 
from IReines. Johnson, fc GossI ( 2008l l and assume the therm al emission follows oc and 

non-thermal emission as oc (e.g. iBaars et al.Nl977l ). The thermal emission is, however, 

dependent on both the size and density of the emitting region as S'jy^thermai oc IkTu^Tyj(?, and 
the models are thus under-constrained. To better constrain the fit, we can adopt the approxi¬ 
mate size of the 3.6 cm radio continuum emission region. The 3cj contour suggests a radius of 
^ ^50 pc - however a smaller siz e would be derived if only the unresolved peak emission was used 
(jReines. Johnson. &: Gossll2008l . also seen in Figured]). 


We find the model with the lowest by both setting the radius and leaving it as a free 
parameter, as shown in Figure O By imposing a radius of the adopted value of 50 pc, we find the 
density is rie = 21.5 cm“^ and the non-thermal contributions at 1.3 cm are 14% ± 15%, the band 
least contaminated by non-thermal emission. If the radius is unconstrained, we find the thermal 
radio emission results from a region of radius of r= 2.3 pc, a density of = 2.0x10^ cm“^, and 
non-thermal contributions at 1.3 cm of 28% ± 25%. The fit is under-constrained for both models 
iXred 0.5 and ~0.1 for a set and free radius) and the equality of the fits cannot discriminate 
between the two conditions and rather express both as possibilities. 


Thus, the radio data is consistent with both a large, low-density HII region typical of Giant HII 
regions (GHRs) and a very dense, compact region typical of UDHHs (or, more likely, a combination 
of the two, although the data do not permit a more sophisticated model). S26 clearly has associated 
extended thermal radio emission, thus there may be a low filling factor of small dense regions or a 
single very dense region within a larger, less dense S26 HII region. This type of scenario has been 
observed in Mrk 996, where the extremely dense (ne~ 10® cm ^) nucleus is sur rounded by a less 
dense (~10^ cm“^) star forming region, along with WR stars ( Telles et al.ll2014l '). We find the flux 
at 1.3 cm of Sj/q.a cm~1360 /rJy is roughly representative of the thermal radio emission. While S26 
is not solely a thermal source, decomposition shows the flux at 1.3 cm can be considered primarily 
thermal. 
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While the majority of the thermal radio emission may be emitted by an HII region surrounding 
the cluster, the Wolf-Rayet populations will also contribute to the observed radio flux. Strong WR 
winds are known to produce thermal free-free radio emission, approximated by the scaling relation 


S'v oc 


M 


(jCrowtheiil2007l ) . Using this relation, we can estimate the expected radio flux density contribution 
from the Wolf-Rayet stars in S26 to identify whether the HII region or the WR winds are the likely 
dominant emitter. Accounting for metallicity, we can estimate the expected radio flux density by 
scalin g from galactic WR stars observed with the Very Large Array (ICappa. Goss, &: van der Hucht 
20041 ). Exactly how mass-loss rates depend on environment is not well known, although mass-loss 
rates at low metalliciti es are observed to be lower than higher metallicity environments as M-wn oc 

2QQ^. Clumping in the winds can additionally reduce 


70.7 


and Mwc oc (ICrowther fc Hadfield 


mass loss rates by factors of two (jCrowtherl 1200711. Taking a generic WC star for example, WR 5 
(WC 6 star) with = 0.20 mJy at 3.6 cm (|Cappa. Goss, k, van der HuchtI 120041 ') would produce 
S 3.6 cm~5.2 X 10 ® mJy at the distance and metallicity of S26 in NGC 4449. Similarly, the radio 
flux density of a WN star could produce 83,6 cm~5.6 x 10“’^ mJy (e.g. the WN 7 star WR 100). 
Summing the total population of WR stars in S26 and scaling to 1.3 cm, we might expect a flux 
density of around Si/, 1.3 cm 9 xl0“® mJy due to the WR winds, which compared to the observed 
flux of Si/_i .3 cm~l-360 mJy, is negligible. Therefore, the thermal radio emission at 1.3 cm resulting 
from the WR winds is 1:10“® to the thermal radio emission observed in S26 and the majority must 
be resulting from free-free emission from the HII region ionized by the stars. 


6. Dissecting the Dust 

We utilize archival data of the galaxy NGC 4449 from near-infrared to far-infrared wavelengths 
to construct the spectral energy distribution (SED) of the massive star cluster S26. The infrared 
emission from S26 is clearly one of the dominating sources in its host galaxy, shown in Figure O 
and at the longest wavelength (500 ^m) is comparable in brightness to the nuclear emission. 


6.1. Photometry 

Embedding material surrounding S26 is evident throughout the rich archival dataset of NGC 
4449 and can be characterized through photometry and the construction of an SED. An SED pro¬ 
vides key information on the heated dust properties and is crucial in describing the evolutionary 
phase. However, the wavelength coverage from a few to hundreds of ^m, along with source com¬ 
plexity as well as corresponding worsening resolution, necessitates careful treatment of the different 
datasets to preserve the shape of the SED. 























- 15 - 




Fig. 5.— Identifying th e ther mal radio component in the radio SED of S26, as observed by 
Reines. Johnson. &: GossI (|2008l h with a fit to the observed fluxes of the expected emission 


assuming thermal emission goes as oc v~^'^ while non-thermal emission as F^, oc . Left: 
Constraining the fit such that the radius is set to 50 pc, as observed by the irregular 3a contour 
of the radio continuum, we find that the flux has non-thermal contributions of 14% at 1.3 cm. 
Right: Allowing the radius to be a free parameter (although under-constraining the fit) results 
in an increased non-thermal contribution of 28% at 1.3 cm from a region of radius 2.3 pc. Color 
version available online. 


We perform photometry using two apertures to evaluate both the resolved emission and the 
total emission, which additionally will provide insight into any radial trends; the apertures are 
overlaid in Figure [Hand explained below. Apert ure photometry is performed using the IDL proce¬ 
dure SURPHOT ([Reines. Johnson, fc Gos3l2nn8l ). The uncertainties are dominated by background 
subtraction, and these are estimated empirically by calculating the standard deviation of the fluxes 
measured using different backgrounds. 


Using native PSFs, we use a small extraction aperture for images where S26 is resolved from 
nearby sources (images centered at 100 /rm or less). We adopt a small circular aperture with a radius 
of 5” (~100 pc; see Figure[3]), which sufficiently excludes nearby sources while including the resolved 
emission from S26. We follow the IRAC Instrument Handbook for Calibration of Extended Sources 
http://irsa.ipac.calt ech.edu/data/SPITZER/ docs/irac/iracinstrument handbook/29/#_T oc296497401 
as well as the procedure in Remv-R.uver et al. ( 2013 ) and the encircled energies in lBalog et al. ( 2014 1 
for extracting fluxes from PACS data, in agreement with http: //herschel. esac. esa. int/Docs/HerschelUG/HUG 


The total emission is measured with a large aperture over the complete dataset. However, 
resolution matching is necessary due to the large change in PSF with increasing wavelength for 
the IR datasets. Therefore, we construct the total emission SED using a large aperture extraction 
region (a 15”~300 pc radius) on images convolved to a common resolution (SPIRE 500 iim with a 
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Fig. 6.— An SED of the optical to far-IR emission from the massive star clust er S26 with dust ari d 
grain models over plotted for comparison. Optical spectra (green line) is from iReines et al.l (1201(1 ). 
and a Starburst99 model (dotted green line) clearly demonstrates an abundance of IR emission to 
the expected stellar emission. The IR photometry and dust and grain model fitting utilize two 
apertures, explained in further detail in Section [6) A small 5” aperture (blue squares used to 
extract emission from images with native PSFs, and a large 15” aperture (black squares) used to 
extract emission f rom images convolv ed to the resolution of SPIRE 500 /rm. The best fit dust and 
grain model from iDraine &: Lil (|2007l l are plotted for the corresponding aperture (large or small), 
where the dotted line shows the model spectra and the empty circular dots show the photometric 
fluxes. The SED clearly shows strong PAH features and large amounts of dust surrounding S26. 
The change in PAH emission in the IRAC bands at 5.8 /rm and 8.0 /xm and the shift in the peak 
IR emission can be seen between the two IR aperture SEDs, plausibly due to increased exciting 
radiation and destruction of PAHs towards the cluster center. 


PSF of 36”) using convolution kernels provided bv lAniano et al.l (|201ll ). The aperture correction 
for images convolved to the SPIRE 500 /xm PSF and extracted with the large aperture are found 
by applying the same photometry procedure to the corresponding PSF and the same extraction 
apertures, as the region size is sufficiently small enough to approximate a point source at this PSF. 


6.2. Spectral Energy Distribution of S26 

The spectral energy distribution (SED) of S26 is shown in Figure [U the shape and peak result 
from large amounts of heated dust. For comparison, we plot the expected stellar emission from 
the Starburst99 model (see Section 17.1. Ij) of a cluster of similar mass undergoing an instantaneous 
starburst. It is clear that the infrared emission cannot be produced solely by stellar sources. Instead, 
this emission is produced from a combination of heated dust partially embedding the cluster and 
Polycyclic Aromatic Hydrocarbon grain (PAHs) emission in the photodissociation region (PDR) 
surrounding the ionizing stars of the cluster. Strong evidence for PAH emission is observed as an 
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excess from the IRAC 5.8 /rm and 8.0 /rm - viewed as a bump between the IRAC 4.5 /rm and 
MIPS 24 /rm points. The infrared emission is extracted using two different methods due to the 
drastic change in resolution at increasing wavelengths, as explained above. Not surprisingly, the 
SED extracted from the large aperture results in more emission than the SED resulting from the 
smaller aperture. 


6.3. The Best Fit Dust and Grain Model 


SEDs can be quantitatively analyzed to describe the emitting material: to characterize the 
dust, PAH components, and exciting radiation, we compare d ust and grain mode ls of a gas and 
dust mixture heated by a distribution of starlight intensities from iDraine &: Lil (j2007l ) to the infrared 
SED of S26. The starlight intensity is given by a scaling factor U such that the energ y density per 
unit fr equency is Uy = where the interstellar radiation held is from 


Mathis et al. 


(Il983l i. The fraction of the dust mass that is exposed to a distribution of intensities from Umin to 
Umax is described by the parameter 7 ; therefore (I- 7 ) is the fraction exposed to Umm - The fraction 
of the total dust mass that is in PAH particles is qpah and expected to be low at the low metallicity 
of NGC 4449. 

The best ht models to the S26 SEDs are found via a test between the observed photometric 
data for S26 and a model’s photometric hux at corresponding bands. The characteristics of the 
best ht models are listed in Tableland the hts are shown in Figure [H The best ht models suggest 
that the exciting radiation in S26 is strong, with a maximum starlight intensity of Umax =10®, 
which is similar to starbursting galaxies. The minimum exciting intensity may increase towards the 
cluster center, as Umin =12.0 for the best ht model to the large aperture SED and Umin =25.0 for 
the best ht m odel to the small apertu re SED. Compared to a global estimate of 2% PAH emission 
in NGC 4449 (jKarczewski et al.ll2013l ). the best ht models may additionally indicate this changing 
environment towards the center of S26. The large aperture SED (extracted from a radius of ~300 
pc) best matches a model with high PAH emission (3.19%) and nearer to the center of the cluster, 
the small aperture SED (extracted from a radius of ~100 pc) best matches a model where the 
PAH emission is reduced to 1.12%. Thus, the dust grains may be being destroyed due to a higher 
radiation held. 


6.4. Dust Mass and Star Formation Efficiency 

Dust is an integral component to the natal cluster environment. Understanding the dust 
mass that is surrounding a region is important in determining the evolutionary stage, providing 
comparisons of similar regions to formulate an observational picture, and for extrapolating to 
estimate the gas mass. The ratio of the stellar to gas mass serves an indicator for the degree of 
removal of natal material. Utilizing results of the dust model fitting and the observed photometric 
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Table 3. SED of S26: Photometry and Model Parameters 



Small Aperture 

Large Aperture 

Band 

Flux Density (mjy) 

Flux Density (mJy) 

IRAC 3.6 fim 

1.98 (0.21) 

9.1 (3.1) 

IRAC 4.5 fim 

2.07 (0.18) 

8.8 (3.4) 

IRAC 5.8 fim 

6.46 (0.42) 

45 (16) 

IRAC 8.0 fim 

17.24 (0.93) 

139 (52) 

MIPS 24 fim 

174 (17) 

355 (88) 

PACS 70 fim 

1180 (180) 

3760 (770) 

PACS 100 ^m 

1090 (150) 

6100 p600) 

PACS 170 Atm 


4690 (990) 

SPIRE 250 fim 


1640 (250) 

SPIRE 350 /rm 


720 (120) 

SPIRE 500 /rm 


290 (30) 



25.0 

12.0 

U max 

10® 

10® 

7 

0.08 

0.03 

Qpah 

1.12% 

3.19% 


^Notes: The upper part of the table presents the flux densities 
extracted from the corresponding aperture. The lower part of the 
table presents the param eters determined by the best fit models 
from lPraine Sz Lu {2003), which describe a gas and dust mixture 
heated by a distribution of starlight intensities, U, where the local 
starlight intensity is U=l. The fraction of the dust mass exposed 
to Umin is I- 7 . The emission due to the PAH particles is given 
as QPAH- 
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fluxes, the dust m ass surrounding the cluster S26 can be estimated. 


Following the method of 
Draine &: Lil ( 20071 1 . we And a total dust mass of 2.3 xlO® Mq as measured with the large aperture 
(with uncertainties of 15% due to flux uncertainties alone). Due to the resolution at the long 
wavelengths necessary for the dust mass estimate, we only have a measurement over this spatial 
scale, which is ~300 pc. 

This total dus t mass, which is sam pled over a large region, is roughly three times the stellar 
mass estimated bv iReines et al.l (2010|) of 6.5x10^ Mq, which corresponds to a region of r=3”.3. 
However, we can scale the stellar mass to the same large 15” aperture used for the dust mass 
estimate. We determine this scalefactor by comparing the flux measured in an archival Hubble 
Space Telesc o pe I-b and image extracted with the r=15” circle (as above) to the value published in 
Heines et al.l ( 20111 ). We And the scaled stellar mass in S26 over the large aperture is 3.0 xlO^ Mq. 
To estimate the total gas mass surr ounding the c luster as well, we assume a dust-to-gas ratio 
suggested by the best fit model from iDraine &: Lil (j2007l i of 1/130, which gives a gas mass of 3.0 
X 10^ Mq. The uncertainties are at least 50% based on the range of acceptable models a.nd ad opting 
a dust-to-gas ratio. For comparison to the host galaxy, modeling bv iKarczewski et al.l (j2013l ) found 
a global value of dust-to-gas ratio equal to 1/190 in NGC 4449. Regions of higher star formation 
events such as in LIRGs s how average value s of the gas to dust ratios of 120 ± 28 (Luminous 
Infrared Galaxies - LIRGS: IWilson et al.ll2008l ). in excellent agreement with the model value found 
for S26. By adopting the above total gas mass, the star formation efficiency (SFE) of the entire 
S26 region can be estimated. 

A gas mass this large implies a low global SFE of 1% for the region of r^300 pc surrounding 


Krouoa et al. 

2001 

). However, low SFEs can be observed, such as 5-10% over entire molecular 

clouds ( 

Williams &: McKee 

1997 

), and the size scale over which SFEs are measured can have a 


major role. Additionally, the extremely low inferred SFE of S26 may result from several factors 
that we cannot discriminate between. The extraction region for the large aperture SED is roughly 
35 times the area used to estimate the stellar mass of the cluster, and thus may be including dust 
and gas beyond the S26 structure. Alternatively, we cannot rule out if S26 is done forming stars 
or if the massive stars are disrupting or inhibiting further star formation that might be occurring. 


7. Discussion 

It is evident that the massive star cluster S26 in NGG 4449 is emerging from its natal cocoon, 
as seen through multi-wavelength evidence: 1. thermal radio emission from ionized gas indicative 
of a young massive star cluster, 2. unique optical features produced by the evolved massive Wolf- 
Rayet stars producing large amounts of mechanical stellar feedback, and 3. strong infrared and 
PAH emission from a PDR and heated dust surrounding the cluster. Here we investigate physical 
conditions pertaining to the massive stars within the cluster, which may be driving or contributing 
to this cluster’s emergence from its birth material. 
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log(Time[yr]) 


Fig. 7.— The estimated age versus the observed WR and 0-star populations, and WR/0 pop- 
ulation ratio, of the cluster S26 (data points). The evolutionary synthesis model of Starburst99 
(jLeithereiill999f l. which utilizes a Kroupa IMF with upper and lower limits of 0.1 - 120 Mq 
plotted (dotted lines) and further discussed in 17.1.11 Color version available online. 


IS 


7.1. Massive Star Populations of S26 in Context 
7.1.1. Estimates through Starburst99 


Models and population comparisons may show evolution timescales, metallicity or environ¬ 
mental effects, or deviations from the expected IMF, and thus indicat e processes that are affecting 
the star formation. Using Starburst99 v7.0.0 (jLeitherer et al.lll999l) . the starburst properties of 
S26 are modeled by adopting an input metallicity of Z = 0.008 (see Section 1,1.2h and simulating 
an instantaneous burst of star fo rmation. The Geneva evolutionar y tracks with high mass loss and 
Pauldrach/Hillier atmospheres ( Smith. Norris, fc Crowther 2 OO 2 I 1 are used with a Kroupa IMF. 


The Starburst99 models are used to determine the age (Section 13.411 and provide an evolutionary 
comparison for theoretical massive star populations, as shown in Figure [71 To produce comparable 
massive star popu lations to the observed data, we scale the initial cluster mass to 6.5 xlO'^ Mq 
([Reines et al.ll201Cll l and use an upper mass limit in the IMF of 0.1 - 120 Mq. The observed WR/0 


population ratio in the massive star cluster S26 is consistent with these predictions, but is slightly 
higher than expected, as shown in Figure [71 If the age were a little bit larger, this discrepancy 
would go away. In general, we expect fewer Wolf-Rayet stars at low metallicities like S26. 
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12 + log{0/H) 



Meynet & Maeder (2005) 

- Geneva model (with rotation) 

- Geneva model (no rotation) 

• Galactic values 

• Cluster/complex values 


Fig. 8.— A compilation of observed population ratios in nearby galaxies compared to differ¬ 
ent predictions. Star points show the observed ratios in the m assive star cluster S26 in NGC 
4449. Lines show model predictions of (jMevnet &: Maedeii l2005l l (dotted; including stellar rota¬ 
tion) and computed results for new models from the Geneva group (solid) for z=0.006 (Neugent 
et al. 2012) and z=0.014 (Georgy et al. 2012) in blue and green. The symbols roughly in¬ 
dicate the method of determining the populations: squares show individually, spectroscopically 
resolved populations and circles show ratios similarly measured to th is work (open are correct ed for 
completeness). Golor version available on line. Refe rences: 1-IC 10 (IMassev &: HolmesIl2002lL 2.3- 


IG 4 662 A1,A2 ( Crowther fc Bibbvl l2009lL 4-LMG(Neugent et al.N2012l ). 5-M 31 (jNeugent et al 


2OI2IL 6.7.8-M 33 inner, midd le, outer ( Neugent et al.l I2OI2IL 9/9a/9 b-M 83/M 83-74/M83 


Hadfield fc Growthei 


^ - 1 ■ 1 - - - 7 ~ -7 T p —■ ■ = ■' " I I - ' - ' - 

Hadfield fc Growthei 2008). 10-Mrk 996 nuclear ( Telles et ah 2014), ll-MW ( Georgy et al 
201^1. 12-NGG 300 (Ribbv fc Crowtherl boiol l. 13-N GG 1140 (Imq]] et abl boOTtl. 1 4-NGG 1313 


in WR/0 : red filled circle), 16-NGC 3049 ( Gonzalez 


(Sidoli et al 


20071 ). 15-NGG 1569 SSGA ( Gonzalez Delga do et al.1 119971 ) (under S26 


A1,A2,B (Hadfield fc CrowtherlI2OO6), 20-NGC 5068 (Bibby &: Crowtherl I2OI2I ). 21-NGG 5253 B 


( Bibbv &: Growtheiil20ld L 24-SMG ( Neugent et al.ll2012lL 25-Tol 89 ( Sidoli et al.ll2006lL 


Delgado et al 


200: 


2004lL 22-N GG 6822 (IMassev fc Johnsonl 1 199811 . 23/23 a-NGG 7793/NGC 7793 R34 


ik 


17,18,19-NGC 3125 


7.1.3. Wolf-Rayet Population Trends with Metallicity 


It is necessary to additionally compare the massive star populations of S26 to observations 
of other regions and models that are in different environments. The number of Wolf-Rayet stars, 
specihcally comparisons between the subtype of the WRs, is heavily dependent on the metallicity. 
The ratios of Wolf-Rayet populations WC/WN and WR/0 stars in the local group are observed 


to roughly decrease with host galaxy metallicity (e.g. lMassevlll996l ). Figure [8] shows evolutionary 


predictions against observational data of nearby gala xies and S26 in NGC 4449 . Evolutionary model 
predictions of the Geneva group are sho wn, includinglMevnet &: Maedeii (|2005l l , for initially rotating 
single stars, and new results presented in iNeugent et al.l (1201211 for Z=0.006 and lGeorgy et al.l (j2012l l 
for Z=0.014, with and without including stellar rotation. We convert from theory predictions at 
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each metallicity Z to 12+log(0/H) assuming a simple scaling relation and 12+log(0/H)Q = 8.69 
( Asplund et al. 2009 1. 


Similar plots display i ng m etallicity trends of Wolf-Rayet population s are shown by the survey s 
fror niBibbv fc Crowtherl (120101 1 as well as a very thorough sampling by iNeugent fc MassevI (120 111 ) 
and iNeugent et al.l (120121 ). Here, we compare the inferred stellar populations in S26 to a larger 
sample of other massive clusters and galaxies found in the literature (references given in the caption 
of Figure [8|). In several cases, observational datasets have been evaluated for completeness and/or 
individual stars have bee n spectroscopically confirmed as Wolf-Rayet stars fe.g. INeugent et al.ll2012l : 
Bibbv fc Crowther 2010 1. 


The inferred subtype WC/WN ratio of S26 roughly agrees with predictions (Fig. [8]) in the low 
metallicity cases, especially in the cases that have been corrected for completeness. The observed 
WC/WN ratio in S26 in NGC 4449 does not appear unusual in comparison to other galaxies. 

In contrast to the relatively well-behaved WC/WN ratio, scatter is seen between the WR/0 
ratio observations and predictions at high WR/0 ratios that has received attention in only a few 
cases. This is rather surprising, as recent comparisons of the observed WR/0 ratio to predictions in 
the so lar neighborhood w ere used to suggest that single star evolution may only account for 60% of 
WRs (jOeorgy et al.ll2ni2l ). As Figure [8] shows, a few observed regions populate an area on the plot 
of low metallicity yet high WR/0 ratios. These values often correspond to SSCs rather than values 
from integrated regions across galaxies (although the distinction becomes unclear for the smallest 
galaxies). For instance, points plotted for NGC 3125 and Tol 89 are known intense star-forming 
regions. There was debate over how high the intensity of the starburst and WR populations in NGC 
3125 really are ( Schaerer. Contini. fc Pindao Il999 ; Chandar et al. 2004 ; Hadfield &: Crowther 2006 ; 
Wofford et al.ll2014l h and in this case the high WR/0 ratio was broadly discussed. Tol 89 is one of 
the brightest known GHRs and has bee n resolved into individual compact clusters, four of which 
contain WR clusters (jSidoli et al.ll2006l b These regions clearly display large WR/0 ratios which 
might then suggest some sort of extreme star formation indicator. Most recently, the unusually 
high WR/0 r atio in NGC 3125-A l has be interpreted to suggest the upper mass limit of the IMF 
is > 120 Mq ( Wofford et al. 2014 1. While the massive star cluster S26 in NGC 4449 is observed in 
the middle of the observed WR/0 ratios and within the uncertainties is consistent with predictions, 
S26 is currently amid regions of intense starburst nature. 


7.1.3. Additional Considerations 


Many assumptions are necessary to estimate massive star populations, each with caveats that 
could impact the inferred stellar content. First, the 0-star and WR estimates would be incorrect 
if the contamination by Of-stars were underestimated, or if the value of the standard ionizing flux 
Qo f or any of the subtypes were incorrect or changed. In fact, use of the ionizing flux Qo = 


from iMartins et al. (120051 ) results in twice the 0-star population than that of the canonical value 
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of 10^® from iLeitheren ( 19991 ). Additionally, the age of S26 is important in estimating rjo to correct 
from an 07V population to the complete 0-star population and, even more so, for comparing to 


the p redictions of Starburst99. Lastly some WR stars can form from binary systems ([Georgy et al 


2012l i. however we compare to models which include only single star evolution tracks as reliable 
binary tracks are not yet available, which could thus underestimate WR populations at later times. 
Thus, the inferred populations of both WR and 0-stars may be altered with different assumptions, 
which would alter the interpretation. 

While consistent with predictions, Figure [7] shows S26 may have a somewhat high WR/0 
ratio, albeit with large uncertainties. While this discrepancy is not statistically significant, there 
are many physical scenarios that could explain an offset if re al. Firstly, stochastic behav ior becomes 
increasingly important for clusters with masses < 10^ Mq ( Fouesneau fc LanconI 120101 ) and S26 is 
at that limit. A high WR/0 ratio could also be produced with a different IMF or a higher upper 


mass limit than assumed, as in the case of the extreme star formation in NGC 3125 (Wofford et al 


20141 ). Alternatively, the population ratio could be altered by a multiple burst scenario. Lastly, 
if S26 has not yet fully emerged from its natal material, some of the stars may be embedded and 
thus unseen at optical wavelengths. Gomparisons of the ionizing flux seen at optical and radio 
wavelengths show this is unlikel y in S26. The op t ical i onizing flux of S26 is simply estimated 
through empirical relations from ISchaerer Sz Vaccal (199a). Using H/3 at 486lA, the ionizing flux 
is Qo 180x10^® photons s The ionizing flux, as measured by the thermal radio emission, can 
be determined by 


Qlvc > 6.3 X 10 


52, 






0.1 


L 


u, thermal 




GHz lO^^erg Hz 


-1 


([Gondon Ill992l l is QLyc ~ (200±20) x 10^^ s ^ at 1.3 cm. The ionizing fluxes agree within the 
uncertainties, and therefore it is unlikely many 0-stars remain embedded. 


7.2. Impact of the Massive Stars on the Cluster Evolution 

7.2.1. The Potential Ionized Bipolar Outflow in S26 


The morphology of S26 as seen in archival Hubble Space Telescope imaging with narrow band 
filters centered on the ionized lines of Ha, [NHI], and [OHI] suggest that it could be driving an 
ionized outflow on an intermediate scale, dwarfing individual stellar outflows seen in the Milky 
Way yet smaller than a galactic outflow. As shown in Figure El the resolved morphology of the 
central ionized nebular gas in S26 appears bipolar and is evident in all archival HST images taken 
with ionized gas filters, with a size scale of roughly 1” ~ 18 pc. The butterfly morphology of S26 
most resembles that of bipolar HII reg i ons, w hich are a type of UGHH according to the modified 
classification scheme bv iDe Free et al.l (|2005l L As a class, the kinematics of bipolar HII regions 
suggest an ionized outflow directed by a central source. S26 appears the most morphologically 
similar to the Galactic bipolar HII region S106, which is surprising given that S106 spans roughly 
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Fig. 9.— An HST Ha image zoomed in on S26, showing the morphology of a bipolar ionized 
outflow. The green line shows a 3 cj contour from 3.6 cm continuum, as in Figured) 


0.006-0.009 pc across and is ionized by a single 0-star (jChurchwelll l2002l l . quite different from 
the massive stellar populations contributing to the nebula in S26 in NGC 4449. Although vastly 
different spatial scales, both S26 and S106 display a dark lane that bisects the hourglass nebula. 
In S106, the dark lane is thought to result from a combination of shadowing and pro tection of an 


inner disk and high column-density, warm gas on the edge of the molecular cloud (I Simon et al 


20121 ). However, many factors may contribute to this apparent feature in S26 given its distance 
and complicated environment. 

The intermediate size of possible outflow driven by S26 is quite intriguing if the outflow is 
confirmed. In the Milky Way, ionized outflows from protostars are commonly on the order of a 


parsec in size (jBachiller 1 19961 ). Alternatively, SSCs may contribute galactic winds or outflows 


that typically span several kpc yet do not appear as localized as S26. The starburst in NGC 1569, 
including several SSCs, appears to be driving a (uncollimated) massive outflow seen as diffuse X-ray 
spurs corresponding t o well-known Ha hla, ments with a high-velocity expanding component that 
is over 2 kpc in size (iHeckman et al.lll995l ). Morphologically, S26 is more similar to the bipola r 
superwind of M82, possibly driven by dense clustering of SSCs (e.g. IWestmoauette et al.N2007l b 
or perhaps the superbubble off the nucleus of N GC 3079 that is 1 kpc in diameter and powered 
either by an AGN or starburst (jCecil et al.ll200lh . These outflows are clearly much larger than the 
ionized region in S26 - thus If S26 is a verified bipolar outflow, it may provide insight into potential 
precursors to, or alternatively failed, galactic winds. 
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7. 2.2. The Importance of Winds from Evolved Stars 


For some clusters like S26, it may be possible that the mechanical luminosity due to strong 
stellar winds of an evolved massive star population is the tipping point in the cluster evacuation 
process. In dense, high-pressure environments, the combination of typical feedback processes from 
forming stars (such as radiation pressure) may be insufficient to clear a cluster before a supernova. 
Yet massive stars, on the main sequence and more evolved, could produce the necessary additional 
mechanical feedback via wind, in combination with the other forms of feedback, to completely 
clear a cluster prior to a supernova explosion. However to occur before a supernova event, the 
removal would have to occur over a very short timescale. It is conceivable the enhanced wind 
phase of evolved stars could drive this evolutionary transition even more efficiently. Thus, perhaps 
in certain environments, the question is whether the wind mechanical feedback is required, rather 
than the dominant feedback mechanism, to clear the cluster before a supernova, which will certainly 
alter the cluster’s ability to survive. 


Massive stars driving cluster evolution via winds have been observed; [Gilbert fc Grahaml (120071 ) 
suggest that ELCs are evacuating their surroundings through wind, finding evidence for outflows 
as the HII gas is not bound. Yet, the extent to which the mechanical luminosity due to stellar 
winds, particularly those from evolved stars, may contribute has not been fully investigated in the 
SSC regime, yet it must be substantial. The interaction between the photo-ionized HII region and 
the stellar wi nd bubble strongl y affect the morphological evo lution even in the case of a single 
massive star (|Frever et al.l 120031 1. ISilich &: Tenorio-Taglel (j2013l l show that radiation pressure on a 
wind-driven shell of a cluster becomes negligible after 3 Myr and highlight the importance of the 
mechanical feedback. 


The first step in answering whether the winds from evolved stars are an essential component in 
the feedback scheme driving cluster evolution is to investigate if the massive star habitants of S26 
could fully remove its natal material solely via winds. We compare the binding energy of the initial 
cluster in S26 and to calculations of the energy input by the massive star population, highlighting 
the additional boost given by the evolved WR stars. 


We first approximate the binding energy of the cluster as a rough estimate of the energy 
necessary to remove embedding dust and gas. The region corresponding to the 3.6 cm radio 
continuum in S26, which includes the nebular emission with the bipolar outflow morp hology, is the 


regio n being cleared out, called the ‘core’ here. The stellar mass of the core is known (jReines et al 
20101 . 6.5x10^ Mq; ); yet because the dust and gas mass are estimated for S26 over a much larger 


region, we consider a total core mass produced by a range of SFEs. We approximate the total mass 
in this central region by including the gas and dust as twice the stellar mass (a roughly 50% SEE) 
and, closer to what is observed, using a SEE of 5% (jWilliams fc McKeelll997l ). We find the binding 
energy of S26 is Ebmd = 8.8 x 10^® erg (lOpc/r) to 8.8 x 10®^ erg (lOpc/r) respectively, where 
the radius is normalized to a canonical value of 10 pc. Pressure contributions from the outside 
ISM beyond S26 contribute a negligible energy threshold to overcome when expanding the cluster 
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- however the large surrounding dust and gas mass have not been considered. 


We compare the binding energy of the cluster core to the effective energy from the cumulative 
mechanical luminosity output by the winds from massive stars as L^md ~0.5 Mv'^. We assume a 
conservative 1% efficiency in the transfer of the mechanical lumi nosity into the surrou nding material, 
interpreted over integrated 0-star lifetimes from simulations of iFrever et al.l (120061 ). As S26 is ~ 3 
Myr old and thus the WR phase has likely just begun, we adopt a conservative timeframe for the 
WR phase of S26 to be the average WR ph ase of 0.3 Myr for an individual star. We adopt typical O- 
star c haracteristics of M = 10“® Mq yr“^ ( Massev 2003) and Voo = 2000 km s“^ ( Kudritzki fc Pu1^ 


2000|) for an 07V star (along with r]o to estimate the general 0-star population, as in Section [12]) • 


We a .dopt typical WR star mass-loss rates and terminal velocity values of WN7 and WC4 stars 
from ICrowtheii (j2007l ) , and treat any metallicity effects as in previous sections for the WRs and as 
Mo oc for the 0-stars ( Mokiem et al. 2007). 


To first order, assuming a single starburst, S26 may have been completely cleared by the O- 
star winds alone - although the fact that the cluster is still partially embedded indicates this has 
not yet happened or suggests the SFE in the core is indeed low. Over 3.1 Myr, the estimated 
cluster population of 237 0-stars would contribute 8.8x10^® erg through winds. However, 18 of 
these massive stars have evolved into WR stars, which contribute through much stronger winds. If 
we account for increased mechanical luminosity for the number of inferred WRs over the average 
0.3 Myr long WR phase of an individual star, the total massive star population would output 
9.9x10^® erg over the same 3.1 Myr - larger than the binding energy for an initial cluster with a 
radius of 10 pc with a high SFE. Thus, the massive star population has likely contributed enough 
mechanical luminosity alone to have fully cleared the cluster, especially when WR contributions are 
considered. S26 likely has not done so because of the surrounding material or because the SEE is 
low. We hypothesize the pre-supernova mechanical luminosity is crucial in the evacuation process 
in this case, especially as the removal is not instantaneous and thus more likely to ensure cluster 
survival. 


7.2.3. Describing the Evolutionary Phase of S26 


We have examined S26 to investigate this seemingly short-lived yet critical stage in mas¬ 
sive/super star cluster evolution, and we now identify its place among the major classes of HII 
regions. In Eigure m we plot the observed characteris tics of S26 on the size-density relation of 
extragalactic HII regions from [Hunt &: Hirashital (|2009j) . The optically derived density, adopting 
a size of 50 pc as suggested by the 3cr contour of the 3.6 cm radio continuum that contains the 
optical cluster, puts S26 among the ETC regions. The two scenarios consistent with the radio data, 
described in Section (5] show that S26 exhibits characteristics similar to both GHRs or UDHHs, 
depending on whether the radius is allowed as a fr ee parameter or set at 50 pc. In the context of 
the evolutionary classification scheme described bv lWhitmore et al.l (j20I4l ). S26 lies on the border 
between Stage 3 (emerging cluster) and Stage 4 (young cluster). 
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Fig. 10.— The optic al and radio estimated properties of S26 plotted on the extragalactic HII region 
size-density relation (jGilbert &: GrahamI 120071 : iHunt &: Hirashitall2009lh As discussed in Section O 
the observed radio emission can be modeled with an input radius of 50 pc or with the radius as a 
free parameter. Golor version available online. 


7.3. Similarities to 30 Doradus in the LMC 


We compare S26 to 30 Doradus in the LMG located at about 50 kpc ( 

Walker 2012'). whose 

dense stellar core R136 is the ‘orototvoe’ core of a SSG (Massev & Hunter 1 

998). In terms of the 

observed stellar populations, S26 more closely resembles the core R136. The S26 and R136 star 

forming regions have almost identical stellar masses (~6xl0^ Hunter et a 

. (1995) for R13fi and 


fti6.5xl0^ Mi^ lReines et ahl ([20101) for S26). S26 hosts abou t two times the mas sive star population: 
R136 contains ~121 massive stars within a 4.7 pc region ([Hunter et al.lll995l i and S26 hosts 237 


massive stars in a region roughly 50 pc in size. Although R136 may be more concentrated, a tight 
stellar core of S26 may be unresolved (see Figured]). As discussed in Section[5l the radio data in S26 
are consistent with thermal emission resulting from both a large, low-density region (constraining 
a spectra fit to a radius of 50pc) or a small, dense region (2.3 pc; unconstrained radius), and are 
likely indicative of a low hlling factor of high density regions. 

However, much of the core of 30 Dor would be unresolved if at the distance of NGG 4449, and 
we thus compare the entire systems as well. Within the total 30 Dor region, the VLT-FLAMES 
Tarantula Survey found 722 massive stars in a 150 pc radius (of which 500 have been spectroscop- 
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ically confirmed), resulting in a stellar mass of 1.1 xlO^ Mq ( Doran et al. 20131 ) for 30 Dor itself. 


The scaled stellar mass of S26 is estimated to be 3.0 xlO^ Mq in the large 300 pc aperture, which 
would likewise host ~1000 massive stars following the same s caling. Most im portantly, 30 Dor 
contains a remarkably similar gas mass of (1.3 ± 0.5) x 10^ Mq (IKim et al.ll2003l l over an aperture 
of 200 pc, which roughly matches the gas mass of 3x10^ Mq of S26 with the large aperture of 300 
pc. 

30 Dor proves an interesting comparison to S26 as it is close enough to resolve individual com¬ 
ponents and evaluate the impact of the massive star populations on the surrounding environment. 
Out of the 500 confirmed hot luminous stars in the VLT-FLAMES Tarantula Survey, 31 are WR 
or Of stars that contribute ~50% of the wind luminosity and ~40% of the ionizing luminosity in 


the cluster (jPoran et al.l 120131 ). This region in 30 Dor roughly corresponds to the size sampled by 
the small aperture SED of S26 (~100 pc), and although the optical spectra were obtained for only 
the inner 50 pc of S26, the estimated evolved population of 18 WRs (out of 237 massive stars) is 
similar. Thus, the wind and ionizing luminosity can be extrapolated from the resolved 30 Dor case 
and assumed to be just as significant, if not more so, for S26. 

Comparing S26 and 30 Dor can shed light on how feedback mechanisms may become relevant, 
especially later in evolution. Comparisons of the observed pressure due to stellar radiation, shock- 
heated hot gas, warm ionized gas (HII gas), and dust processed IR radiation indicate that radiation 
press ure dominates with in the central 75 pc region of R136, and HII gas pressure dominates beyond 
that (jLopez et al.ll201lh . Another study that paired observations with photoionization models sug¬ 
gests that the hot X-ray gas instead dominates the mech anics of 30 Dor, and t hat radiation pressure 
does not currently play a major role in the structure (jPellegrini et al.l 120111 1. The importance of 
radiation pressure in at least the initial expansion is clear, however these and many other studies 
do not directly include mechanical wind feedback from massive stars. 

While 30 Dor has undergone a more complicated star formation history and likely undergone 
multiple bursts of star formation over the last 20 Myr (iDe March! et al.ll201ll ). several parameters 
of the 30 Dor region parallel those in S26; for one, the gas reservoirs are roughly the same size, 
stellar feedback is ongoing, and mechanical feedback is largely due to the most massive and evolved 
stars. The structure of 30 Dor is naturally better defined, with the core hosting 121 massive stars 
(~5 Wolf-Rayet), with many more in the vic inity, and is consistent w ith the picture that the most 
massive stars curtail further star formation (jMassev &: Hunteiill998l l. S26 hosts 237 massive stars 
(18 Wolf-Rayet) and appears to be undergoing or ending the hrst star formation event as it is ~ 3 
Myr old. Thus, it seems plausible that S26 is a younger version of 30 Dor. 


8. Summary 

We have presented a detailed analysis of the partially embedded massive star cluster S26 in 
NGC 4449. The main results are as follows: 
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1. We estimate S26 hosts massive star populations of roughly 18 Wolf-Rayet stars and 219 O- 
stars from optical spectral lines by assuming standard WR optical line fluxes for a single WR 
star and a typical ionizing flux for an 0-star. 

2. The massive star population ratio comparing the subtypes WN/WC in S26 is consistent with 
predictions and other observed clusters and galaxies. The population ratio WR/0 is also 
consistent with predictions. If unavoidably large uncertainties are ignored, S26 is among the 
high WR/0 ratios observed for other individual clusters of intense star formation. 

3. Partially embedded by dust, S26 is one of the dominating sources of the IR emission in NGC 
4449. The best fit model to the large aperture SED (300 pc aperture) suggests PAH emission 
that is stronger than the NGC 4449 galactic value. 

4. Dust model fitting to the infrared SED extracted from 100 pc and 300 pc apertures suggest 
the exciting radiation increases and PAH emission decreases towards the cluster center. This 
radial trend suggests the dust grains are being destroyed from within, likely by the massive 
star feedback. 

5. As estimated by the infrared photometry and the best fit model to the large aperture SED, 
the total dust mass of the S26 structure is 2.3x10^ Mq; by assuming a dust-to-gas ratio of 
130, we estimate a gas mass of 3x10”^ Mq. This results in a low star formation efficiency of 
1% over the large 300 pc aperture. 

6. We hypothesize that the mechanical luminosity from the evolved stellar winds in S26 may 
be essential to the emerging process for this cluster, which may ultimately influence future 
cluster survival. Resolved HST images of the ionized gas in S26 display an hourglass nebula, 
which may suggest a possible bipolar ionized outflow. Simple energy calculations suggest that 
the winds from massive stars may be sufficient in clearing the cluster (particularly if the star 
formation efficiency is high in the region being cleared). However, as S26 has not been fully 
evacuated at ~ 3 Myr, the increased feedback contributed by the Wolf-Rayet stars may be 
necessary in clearing out the cluster. 

7. The optical characteristics of S26 match those of Emission Line Glusters on the size-density 
relation of extragalactic HII regions. The radio properties of S26 cannot discriminate between 
the UDHH and GHR extremes, which may imply a low filling factor of high density regions 
within a large, low density HII region. 

8. Similarities between the S26 cluster and the super star cluster analog 30 Doradus in the LMG 
may suggest that S26 is akin to a younger version of 30 Dor. 

It is evident that the massive star cluster S26 in NGG 4449 is undergoing an important evo¬ 
lutionary stage in which stellar feedback is particularly important. We propose S26 may be an 
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example of a short-lived, yet important phase in massive and super star cluster evolution, dur¬ 
ing which the complete evacuation of natal material is aided by the mechanical luminosity of the 
massive stars, particularly evolved Wolf-Rayet stars. 
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